DISCRETE DYNAMICS OF COMPLEX BODIES WITH 
SUBSTRUCTURAL DISSIPATION: VARIATIONAL 
INTEGRATORS AND CONVERGENCE 
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Abstract. For the linearized setting of the dynamics of complex bodies we 
construct variational integrators and prove their convergence by making use of 
BV estimates on the rate fields. We allow for peculiar substructural inertia and 
internal dissipation, all accounted for by a d'Alembert-Lagrange-type principle. 



1. Introduction 

A variational integrator is any recursive rule that allows one to calculate dis- 
crete trajectories from initial data and coincides with the discrete Euler-Lagrange 
equation of some discrete (or discretized) Lagrangian. It is a discrete scheme for 
evaluating (numerically) the dynamics described by some Hamilton or d'Alembert- 
Lagrange principle. A detailed description of the essential properties of variational 
integrators is in [18] (see also [17]). 

In particular, when the potential appearing in the Lagrangian functional can be 
additively subdivided in subsystems, e.g. distinct groups of particles or finite ele- 
ments in space in case of continuous problems, it is possible to select asynchronous 
time discretizations. Asynchronous variational integrators (AVIs) proposed in [9] 
and [10] then follow. They preserve the symplectic structure of the original contin- 
uum system. Thus, for non-constant time steps the energy evaluated numerically 
oscillates around its average value. 

Although problems in evaluating individual trajectories are the same as for other 
traditional alghoritms, AVIs perform better in determining time-averaged quanti- 
ties. Artificial numerical damping is not introduced. For a desired error value the 
computational cost is lesser than the one of other methods. AVIs display analogies 
with subciclye methods (see [2] and [4]) and other methods in molecular dynamics 
(appropriate comparisons are in [10]). 

AVIs have been introduced in [9] for the dynamics of elastic simple bodies. Their 
convergence in time has been analyzed variously: A repeated use of Gronwall in- 
equality on the discrete Euler-Lagrange equations at a specific node is the basic 
ingredient of the technique introduced in [10] for potentials with uniformly bounded 
second derivative. T-convergence has been exploited for the zero-dimensional oscil- 
lator in [19], some technical hypotheses removed in [TT]. For the linear elastody- 
namics of simple bodies, the convergence has been also proven in [5] along a path 
based on the direct exploitation of the action functional. 
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Here we start the analysis of variational integrators (in the synchronous and 
asynchronous versions) for the linear dynamics of complex bodies, those bodies for 
which the material substructure prominently influence the gross behavior. We con- 
sider cases in which the material substructures display peculiar inertia (additional 
to the macroscopic one) and internal dissipation. The dynamics of such bodies is 
governed by a d'Alembert-Lagrange-type principle. 

We do no treat any specific model of complex materials. We develop our analyses 
in the general model-building framework of the mechanics of complex bodies ([5], 
P"2"]). In this way we do not specify the type of the substructure under analysis. We 
presume only that its essential geometrical features are described by an element of 
an abstract finite dimensional manifold A4, the so-called manifold of substructural 
shapes. We construct first the linearized version of the natural non-linear theory. 
To this aim we need the isometric emebedding of Ai in a linear space. Although 
such embedding always exists, it is not unique, so its choice is a constitutive part 
of the modeling procedure. After establishing the structure of the linear theory, we 
construct appropriate variational integrators and we prove the convergence of the 
discrete approximation in time to the continuous counterpart, the convergence in 
space being assured by standard theorems. 

The technique that we use here is new: it is based essentially on BV estimates 
on the time-rates of the fields involved. 

Convergence in the asynchronous case is proven only for substructures admitting 
quadratic peculiar independent kinetic energy. The expression of the substructural 
kinetic energy is accepted in fully generality in proving convergence of synchro- 
nous variational integrators. Our work prepare the way to a number of numerical 
experiments in various classes of complex bodies. 

Some notations. We use the standard symbols | ■ | for the euclidean norm, 
and both • and (•,•) for the scalar product in some Euclidean space MP, j G N. 
We will not explicitly indicate the dependence on the dimension j, for the sake 
of the conciseness of the notation. The symbol # denotes the cardinality: #T is 
then the cardinality of the set T. By d z we indicate the partial derivative with 
respect to z. We use standard notations for Lebesgue spaces, Sobolev spaces, and 
the space of functions with bounded variation (BV in short). The symbol 
denotes the (d — l)-dimensional Hausdorff measure in M. d . In some computations 
we find inequalities involving constants which depend on the data and on the space 
dimensions, but are always independent of the set of discrete instants. Since it 
is not essential to distinguish from one specific constant to another, we indicate 
all of them by the same letter c, leaving understood that c may change from one 
inequality to another. 

2. Commentary to the dynamics of complex bodies 

In a primitive approach, a body is considered as an abstract set 03 the elements 
of which are called material elements. Each element is considered as the small- 
est piece of matter characterizing the nature of the material constituting the body 
under examination, that is a patch of matter made of entangled molecules or the 
characteristic cell of some atomic lattice. The representation of this set is a model 
of the morphology of the body. It cannot prescind from the selection of the place- 
ments of the body in the ambient space R d , placements that are selected by means 
of one-to-one maps / : 03 — > M d . The generic B :— / (03) represents the 'gross' 
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shape of the body and is assumed to be a bounded domain with boundary dB 
of finite (d — l)-dimensional measure, a boundary where the outward unit normal 
n is defined to within a finite number of corners and edges. When substructural 
complexity arises, description of the structure inside material elements is neces- 
sary. The representation of 95 is then defined by maps x : 25 — > A4 assigning to 
each material element a descriptor of the peculiar features of its inner morphology 
(called for this reason a morphological descriptor), with A4 a finite-dimensional 
cliff crcntiablc manifold. 

The dynamics of a complex body involves then the description of subsequent 
changes in gross and substructural shapes. 

Convenience suggests to select a reference place £>» := /» (05) that can be in 
principle occupied by the body in R d and to obtain from it new places by means of 
sufficiently smooth bijective maps (transplaccmcnts) 

x^y = y(x) = (fof-^ix) eU d , i£8„ 
such that (i) B = y (£>*) is as regular as B* and (ii) the gradient of deformation 

F = Vy £ Horn (T x B*,T y B) 
has positive determinant at x 6 £>*. A time-parametrized family of deformations 

(x,t) ^y = y [x.t) e M d , (x,t) e B* x [0, i] , 

is the motion at gross scale. Sufficient smoothness in time is presumed so that the 
velocity y = (x,t) £ T y i x ^B ~ M. d can be defined as a field over 23*. 
Diffcrentiable maps of the type 

x i — ► v := (xo (x) e M, x e B*, 

define fields of morphological descriptors. Then, a substructural 'motion' is then 
given by 

(x.t) i — > v = v (x,t) e M, (x,t) g £* x [0, t\ . 

Assuming sufficient smoothness in time, the rate of change of the substructural 
morphology in the referential description is v := (x,t). At each x the spatial 
derivative of x i — > v is indicated by N G Horn {T X B*,T V M). 

In a conservative setting the mechanics of a complex body is governed by Hamil- 
ton principle of least action with Lagrangian density 

£ = \p ly| 2 + X (f , v) - pe (x,F, v, N) - pw (y,v) , 

assumed to be diffcrentiable with respect to all its entries, p is the referential mass 
density (conserved during the motion), e the elastic potential and w the potential 
of external actions, all per unit mass. The term \ (u, v) is a real-valued function 
vanishing when v = 0, it is of degree 2 in v and is such that the second derivative 
dl^x exists and is positive definite. Moreover, sup^ (doX " ^ — X) is attained (as a 
maximum) at a unique point in the relevant T v hA. It coincides with the substruc- 
tural kinetic energy when peculiar independent inertia pertains to substructural 
changes (see remarks in [5], [H]). 

The derivatives of C with respect to F and N represent contact interactions, 
namely standard and microscopic stresses, respectively. The derivatives of L with 
respect to y and v are standard and substructural bulk actions, the latter being 
divided additively in external fields acting on the substructure and internal self- 
actions. 
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Dissipative effects may occur. For the sake of simplicity it is assumed here that 
dissipation occurs only within each material element and is accounted for by the 
local self-actions. A concrete example is given by quasicrystals (see [T3], [II])- 

A d'Alembert-Lagrange-type type principle then substitutes the Hamiltonian 
one. One then requires that 

(2.1) 5 CdxAdt+l z° ■ <f> dx A dt = 0, 

JB,x[0,t\ JB,x[0,i\ 

where S indicates the first variation, z v = z v (F, v, v, N) G T*Ai represents the 
dissipative part of the self-actions inside the generic material element and (x. t) \ — > 
4> :— <p(x,t) G T v r Xt fjM is an arbitrary smooth field with compact support. 
To define variations, for s G [0, 1] one considers smooth maps 

£ — *Ve-=y + £<P, £' — >v e , 
with | 6= o (x, t) = 4> (x, t), (x, t) i — > ip :— tp (x, t) a smooth R d — valued map. The 
variation S along (y, v) has then the meaning of the derivative 4^ of C evaluated at 
e = 0. It is also assumed that the self-action z v is intrinsically dissipative, that is 

z° ■ v > 0, 

for any choice of v. Then z° may be of the type 

v • b 

z = r\v , 

with i} the covector associated with v and r\ a positive diffusion coefficient. 

When the fields x i — > P := —pdp£ and x i — > S := p9jve are of class C 1 (£>*) n 
C° (£*), Euler-Lagrange equations are then 

pi) = DivP + b, 

di,x - d v x = DivS — z + (3. 
P is the Piola-Kirchhoff stress tensor, b :— —pd y w the co-vector of body forces, 
S the microstress tensor in referential configuration, z :— pd v e the self-action, 
f3 := —pd v w external body action over the substructure. 

By considering fields over the actual place B, the balance equations above become 

p e a = diva + b e , 



do e X ~ dveX = DivS e - z e + j3 e , 

where 

a := (detFy 1 PF*, b e := (det F)^ 1 b, 
S e := {dctFy 1 SF*, z e :— (det F)' 1 z, j3 e := (det F) _1 (3. 

A requirement of objectivity involving the action of SO (3) on both the ambient 
space R d and M. for the density of the elastic potential e implies that 

skwPF* = e [A*z + (DA)* S) , 

where A = A (v) G Horn (R d , T„M), and e is Ricci's alternating symbol (see [S]). 
Moreover, one may require the invariance of the entire Lagrangian under the action 
of the group of diffeomorphisms of the ambient space and the group of diffeomor- 
phisms of M. into itself. Such invariance requirement opens a path that allows one 
to prove the covariance of the balance laws (see [14]). 
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Boundary conditions can be assigned in terms of places and morphological de- 
scriptors on parts dB* y and dB w of the boundary OB*. Tractions t := Pn and 
zu := Sn, with n the normal to the boundary in all places in which it is de- 
fined, can be in principle assigned on parts dB* t and dB* m of the same bound- 
ary. Mixed boundary conditions require dB* y n <9£>* t = 0, dB* v n 9B* ro = and 
d&Vy U <9B* t = dB w U dB^ = 0. 

A loading device able to prescribe at the boundary the substructural contact 
interactions w is not know in imagination or elsewhere. The sole natural boundary 
conditions in terms of substructural contact interactions is then w = 0. Different 
is the case of Dirichlet data in terms of morphological descriptors. Example is the 
case of quasicrystals in nematic order: surfactants spread along the boundary allow 
to assign a precise orientation to the rod-like molecules in nematic order. 

Initial conditions are given by regular fields (a;, to) — > (yo( x > *o)> u o{ x i to)) an d 
(a;, t ) -> {y (x,t ),v {x,t )). 

3. Linearization 

In the standard theory of elasticity, the space of transplacements can be con- 
sidered as an infinite dimensional manifold. Linearization is developed over the 
tangent space of such a manifold. In such a procedure it is useful to define the 
displacement field 

(x, t) i — ► u := u (x, t) = y (x, t) — x, (x,t) G B* x [0, t] . 

The initial condition is then given by (a;, to) — > Uo(a;,to), and the boundary condi- 
tion is given along dB m = dB^y. 

The condition |Vw| << 1 defines the infinitesimal deformation regime in which, 
essentially, a w P, and the stress may depend linearly on the infinitesimal strain 
e := symVu. 

To linearize the multifield setting described in previous section, an additional 
difficulty arises because M. does not coincide in general with a linear space isomor- 
phic to R k for some k. For this reason, the construction of a Sobolev space of maps 
taking values on M. requires the embedding of the manifold of substructural shapes 
itself in a linear space. We assume that Ai is endowed with a C 1 Riemannian struc- 
ture and the relevant Levi-Civita parallel transport. By Nash theorem an isometric 
embedding of M in a linear space is always available but it is not unique. The 
selection of an embedding becomes then a constitutive ingredient of each special 
model. Once the embedding of Ai has been chosen, one may consider the space C 
of pairs of maps (y, v) can be considered as an infinite dimensional differentiable 
manifold. 

We maintain the assumption of infinitesimal deformation setting, a regime in 
which we can 'confuse' £>* with B, in the sense that if ip indicates a vector field 
tangent to y, then ip w u at any x in £>*. For this reason, from now on we write 
B instead of 23* for the sake of conciseness. Moreover, we also confuse v with 
v e = v o y^ 1 for consequent obvious reasons. 

We consider differentiable fields 

(x,t)< — ► (f (x, t) e T u R d 7 (x, t) G B x [0, t] , 



(z,t)i — ► 4> (x, t) e T U M, (x, i) e B x [0, t] , 
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such that ip and <f> vanish where Dirichlct data are imposed. We continue to write 
v and M as before the embedding. The notation L (A) (u, if) (u, <j>) indicates the 
linearization of A about the pair of maps (u,z>). A superposed bar denotes maps 
calculated at (u,v). At each x and t, the linearizations of the maps x i — ► F and 
v i — ► N about {u,u) are then given by 

L (F) {u,if) =F + Vip 

and 

L(N) (v,<f>) = N + V<f>. 

Previous formulas make sense because of pointwise parallel transport of F over 
curves on R d , and of N over A4. 
The next step is to consider maps 

— >Poc, <r i — >zos, <r i — >5ocr, 

where ? = {F, v, N), and to linearize them about (u,D). The infinitesimal deforma- 
tion setting implies P m a, z w z e , S « S e and we continue to write P, z, S in this 
sense. By presuming sufficient smoothness, we then get 

L (P) (u, ip) (P, <j>) = P + L (P) (V^,V0) , 
L (z) (u, ip) (D, cj>) = z + L (z) (Vip,<j>,V(j)) , 
L (S) (u, (P, 0) = 5 + L< 5 > (Vp,0,V0) , 

where the l/')'s are linear forms of their entries. For the sake of simplicity we may 
assume that (tZ,P) are associated with a (so-called) natural state in which P = 0, 
z = 0, <S = (that is <r = 0, z e = 0, S e = 0), so that 

L(P) (u >¥ >) (P,0) = L( p ) (V^,V<£), 
L(z) (u >¥ >) (p,0) = L« (V^,V<£), 
L(«S) (u >¥ >) (P,0) = L( 5 ) (V^,V0). 

In this case, since the measures of interaction listed above are the partial derivatives 
of e with respect to F, v and N, respectively, the elastic energy is then a quadratic 
form in (Vip,4>,V<p). By some abuse of notation, by taking the common notation of 
standard linear elasticity, we then write 

e = e(U,v,N) 

leaving understood the meaning of the entries, with U = F — Id and Id the identity. 

For the sake of simplicity and without reducing drastically the generality of the 
physical cases covered by our treatment, the assumptions listed below apply. 

(Al) The density of the clastic potential e : R dxd x R fe x R kxd -> [0, +oo) is 
a continuous and coercive quadratic form. There exist constants A, A > 
such that for every (U, v, N) 

\(\U\ 2 + M 2 + |A| 2 ) < e(U, u, N) < A(\U\ 2 + \v\ 2 + \N\ 2 ). 

(A2) The potential of external actions to e C 1 (K d x M fe , [0, +oo)) has quadratic 
growth, i.e. there exists Si > such that for every (u, v) 

|Vto(u,z/)| < 5i(l + \u\ + 
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(A3) The traction t E L 2 {dB, H^ 1 ) is applied on a closed subset dBt of dB. 
Moreover, a homogeneous boundary condition is prescribed for the dis- 
placement u on dB u := dB \ dB t , and for the morphological descriptor v on 
a closed subset dB v of dB. 
Such assumptions cover all special cases of complex bodies that we know in such 
linearized setting (see [TS] for a list of examples) . 

4. Asynchronous Variational Integrators 

The asynchronous variational integrator scheme necessary for the dAlembert- 
Lagrange-type principle (|2.1[) can be constructed easily by following the guidelines 
of [9] and [10]. However, the analysis of the convergence of AVIs in presence of 
an abstract form of the substructural kinetic co-energy is hard. For this reason we 
restrict here the attention to the case in which the kinetic co-energy \ has the form 

where p is a material parameter. The general case is treated later for synchronous 
variational integrators. 

Of course the absence of dependence of \ 011 v itself implies that M. is a linear 
space per se and that an isometric parallelism is selected over A4. The Levi- 
Civita connection resting on a Riemaniann structure mentioned earlier is the case. 
Microcraked soft bodies and incommensurate intergrowth compounds fall within 
this scheme (see [13] and [15] for relevant examples). 

In this case, we then consider the linearized dynamics of a complex body with 
zeroth-order dissipation governed by the dAlembert-Lagrange-type principle 

(4.1) 6A(u,v)[<p,<p]+V(v)[4>]=0, 

which is valid for arbitrary choices of test functions ip, 4>, defined earlier. In (14. 1|) . 
A is defined by 

(4.2) A (B, [t , */] v) := j*' Qf QpH 2 + \m 2 ^j dx-V (B, u, v, t)j dt, 
with 

(4.3) V(B,u,v,t) := / p (e(Vu, u, Vv) + w(u, v)) dx - I t-udH^ 1 
and T> is the dissipation given by 

(4.4) T>(P)[<f>] := r] v-c^dxhdt, 

JBxI 

t] is a positive parameter. 

4.1. Spatial discretization. Discretization in space is obtained by means of stan- 
dard finite elements given by a tessellation T of B, for instance a regular triangu- 
lation chosen to be consistent with the partition of the boundary dB into dBt and 
dB u (see, e.g., 0, [7]). By adopting a triangulation we assume that B is a polyhedral 
domain. The integration nodes, the generic one being indicated by a, are selected 
as the vertexes of the elements K of the triangulation T. 

We consider the space PA (7~) of continuos functions which are linear polynomi- 
als on each K G T. With fixed a bounded open interval I = (to, tf) C R, we are 
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interested in the vector subspace V u of PA (T) ® H 1 (/, M d ) which includes all the 
displacement mappings satisfying u \ qb u = 0. Any map u £ V u has then the form 

(4.5) u(x,t) = Y,M a (x)u a (t), 

where J\f a € PA (T) is the nodal shape function corresponding to node a and u a (t) 
is the value of the displacement at a at time t. Of course we assume u a (t) = if 
a £ dB u . The null condition is not a restriction to our treatment. 

With a slight abuse of notation, we denote by u (t) the vector in ]R D , with D the 
number of degrees of freedom of all nodal placements at time t, relative to the map 
u(-,t). Specifically, D = <i#T, with #T the total number of nodes in T. 

It is possible to check that there exists a constant c-r > 0, depending only on 
the tessellation T, such that 



(4.6) sup \Vu(-,t)\ < c T \u(t)\ . 

When restricted to a generic finite element K, the map u(x, t) in (|4. 5|) is indicated 
by Uk(x, i). Analogously Uxit) indicates a vector in K dx ( <J+1 ) having as components 
the displacements of all nodes of the generic element K. 

We perform an analogous spatial discretization for the field v. In this case 
we have v a (t) £ M. k for every a £ T, t 6 / as a consequence of the embedding 
of the manifold of substructural shapes in M. k . We denote by V„ the subspace 
of PA(T) ® H 1 (i ,lSL k ) which consists of all mappings v satisfying v\qb„ = 0. 
Furthermore, we set M = k#T. 

The spatial discretization At of A is then obtained by restricting to Vr := 
V u x V v the domain of A. More precisely, we consider the semi-discrete action 
At : PA (T) <g> H 1 (l,R dxk ) -> [0, +oo] defined by 



At{u,v) 




if [u, v) £ V T 
otherwise. 



Notice that the nodal pair (u(-),v(-)) belongs to H 1 (I, K- DxM ) for every map 
(u, v) £ Vt- Hence, the expansion in (|4. 5[) allows one to regard At as a func- 
tional defined over H 1 (/,R DxM ). We agree to do that in the sequel in order to 
simplify the notation. 

A further simplification concerns the referential mass density p. In the sequel we 
simply presume that p is constant along the motion. Moreover, we assume that the 
mass matrix associated with the spatial discretization can be expressed in diagonal 
form in some basis. Then, a straightforward computation in that basis leads to 

4r ([to, */];«,!>) := A (K\ {u K ,v K )(t)) dt, 

where, for each K £ T , 

A{(u K ,v K ){%K) := Y, (^K(t)| 2 + ^K(t)| 2 ) -V K {{u K ,v K ){t)), 
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and 

(4.7) V K ((u K) u K )(t)):= I e((yu K) u K) Vu K )(x,t)) dx 



K 

■d-1 



+ / w ((uk, vk){x, t)) dx — I t(x) ■ ujc(x,t) dTL° 

JK JdKndBt 

The diagonalization of the mass matrix leaves invariant the Lagrangian: it can be 
considered, in a certain sense, as a change in observer. 

We call stationary point for At any element (u, v) of H 1 (/, E £)xM ) which for 
any (ip, 0) £ (j, R DxM ) , with ip a = if a e dB u and <f) a = if a e dB v , satisfies 
the d'Alembert-Lagrange-type principle 

(4.8) 6A T (u,v)[<p,<P]+V T (v)[<j>}=0. 

The semi-discretized dissipation XV that appears above is defined by 

where 4> K is the restriction of <fi to K. 

The d'Alembert-Lagrange-type principle (|4.8j) is equivalent to the requirement 
that for any node a £ T \ (<9B U U 9B V ), and any (<p , £ (it a , i/ a ) + Hi (i, R dxk ) , 

if m a := X){K|aeA'} is tlie noda ^ maSS ° f n ° de a and Pa := J2{K\aeK} PK,a, 

we have 

TO a U a (t) ■ <p a (t) dt 

I Y. ( I d F e((Vu K ,UK,Vu K )(x,t)) VAf a (x)dx ) -<p a (t)dt 

Jl {K\a£K} ^ JK ' 

+ Y / d u w((u K ,v K )(x,t))M a {x)dx\ ■ <p a (t)dt 

Jl {K\aeK} KJK J 



( 4 - 9 ) -/ E / t{x)N a {x)dH d - l {x) )-ip a (t)dt, 



and 



jf ■ 0„(<) dt + J Wait) ■ <j> a (t)dt 

= Y, ( d N e((Vu K ,v K ,Vv K )(x,t)) VK{x)dx] -<p a (t)dt 

Jl {K\aeK}^ JK J 

+ [ V (/ d v e{{Vu K ,v K ,Vv K ){x,t))Af a (x)dx) -<t) a (t)dt 

Jl {K\aeK}^ JK J 

(4.10) +/ V (/ d v w{{uK 1 VK){x,t))Na{x)dx) -(t) a {t)dt. 

Jl {K\aeK} KJK J 

Moreover, in case a £ dB u \dB„ only equation (|4. 1 0[) has to be satisfied since u a (t) = 
0, while if a £ 3B v \dB u we have v a (t) = and the sole (|4.9p has to be satisified. 

Results about finite elements for complex bodies imply the theorem below (see 
examples in [15], [16]). 



10 



MATTEO FOCARDI* AND PAOLO MARIA MARIANO 1 ' 



Theorem 1. Assume (A1)-(A3) and consider a family (% n ) of regular triangula- 
tions of B, with m > the mesh size of T m . Let also (u m , v m ) £ Vr m be a stationary 
point for Ax m - The sequence ((u m , v m ))meN converges in H 1 (B x I,M. dxk ) to a 
stationary point of A. 

4.2. Time discretization. For the discretization of the time interval I = (to,tf) 
a partition 9 := {ti} i=0 Ne of [to, tf] with t]\f e = tf is selected. Its size Hq is the 
number 

he ■= max {U+i - ti) . 

i 

Each K s T is endowed with an elemental time set which is an ordered subset 9^- 
of 0. By relabeling the elements we write 



©A = 



1 1 - t K < ... < t K <t R — tfj. 



We assume = UagT0a and, for the sake of simplicity, that ©a R ©a' = {to, tf} 
for any K, K' £ T with K ^ K' . We denote by Tq the maximum of the elemental 
time sizes, namely 

To := maxmaxfti^ 1 — tL). 
k b k 

The circumstance that each finite element can be endowed with a different time set 
is the basic characteristic of AVIs as already mentioned: if one imposes appropriate 
choices of elemental time sets, one may prove conservation of energy in discrete time. 
Take note that in principle such choices could not exist (see related comments in 

For a node a in T, elemental time sets define also the nodal time set a by 

a := |J e K = {t a = t 1 a <...<t^- 1 <t^ = t f }. 

{K\a<EK} 

As a measure of the asynchronicity of 0, we consider the ratio 

max a- e r (maxe K (t 3 ^ 1 - t j K 

(4.11) r e := } 

mincer ( nnn 0f 

and notice that Q 

(4.12) h e < T e < T@h e #T. 

In the sequel we assume that each node a £ T follows a linear trajectory within 
each time interval with end-points that are consecutive instants in a . Such a 
choice characterizes the class of AVIs that we analyze here. Then we denote by 3^e 
the subspace of functions in H 1 (J, E £)><M ) which are continuous and with piecewise 
constant time rates in the intervals in a , and such that u a = if a G dB u and 
v a = if a e dBy. Thus, for each (u, v) £ y&, a <E T, t\ <E Q and t e [t* , we 
set 

Ua " "a (tl) Va (t a +1 ) ~ "a (**) \ 



5k _ 



(?i a ,Z> )(t) 



j.i+1 _ « " j.i+1 _ « 



^Indeed, by definition of time size we have Hq < Tq and tf — io < ?»e#0. For every element 
K we have minx £T {vaine K (t^ 1 ~ ))#8lf < t/ — *0i then by summing on iC it follows that 



min (minft-^ 1 - t 3 „)) < hn#T. 



Inequality l|4.12|l then follows by definition of tq (see <4.11^ ). 
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By following [9], the fully discrete action sum in time is defined for (u, v) £ ye by 

N K -1 

(4.13) At,q{u,v):= E A {K;u k ,v k ), 

KeT j=0 

where for the generic finite element K, A 3 (K; uk, vk) is defined by 

Ai(K;u K) v K ) :=J2 £ ^ ^ " ^ ^ a ^ 
a ^ K {At^ K ^)} 



E E 



" ha \v a (4) I 2 - -^)M(^,^)(4)), 



oe *{i|tie[4>4 +1 )} 
with Vif given by (14. 7p . The choice of the approximation 

{{u K ,v K )(t)) dt « (^ +1 - 

gives rise to explicit integrators of central-difference type and is only one of the 
possible schemes that can be used. 

It is also convenient to define all action sums on the same function space to avoid 
to link the function space itself to the choice of 0. For this reason we extend At,b 
to the whole y := H 1 (l,R DxM ). With a little abuse of notation we set 



At,® (u,u) 



as in P~TU]) if (u, v) e y e 

K +oo a (u,u) ey\y e - 

The discrete stationary points of Ar,e are couples (u, v) g J^e satisfying the fully 
discrete analogous of the variational principle (|4.8p . More precisely, for all test 
functions (</j, (f>) G 3^e H i? 1 (/,K DxM ) there holds 

(4.14) M r ,e(u,iOfa.d+PT,e(*)M = 0, 
where the fully discretized dissipation is given by 

N a 

(4.15) Pr,eWM -EE^)-".^" 1 ))'^" 1 )- 

a£T i=l 

Alternatively, we will refer to the discrete stationary points also as to discrete 
solutions. 

For ^ and 4> as above, by a direct computation we infer 

(4.16) Pr(*)M -EE ^-("-(O " ^(C 1 )) ' 

ogT i=l 

The discretized form of the dissipation has been chosen to be (|4.15p for the sake of 
simplicity. Furthermore, notice that by (|4.15[) and (|4.16[) we have 



\V r {v)[4>]-'DT,e{v 



EE"*^^)-^ 1 )) 



aST i=l 



(4.17) 



N a 



a£T i=l 





< rfTe 


/ v ■ <f) dtdx 




~ 2 


JBxI 
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In Proposition[2]below we will show that given initial conditions (u(to), v(to)) and 
(u(to), i>(to)), the discrete d'Alembert-Lagrange-type variational principle (|4.14[) 
defines inductively a unique trajectory (u, v) G 3^0 ■ 

4.3. Convergence in Time. The following theorem clarifies the stage. 

Theorem 2. Assume (Al)-(AS). Let (Qh)heft be a sequence of time sets for a 
bounded time interval I = (to,tf) such that Tq h — > as h — > +oo. Let also 
(wft(io)j Vh(to)) and (v.h{to),i'h{to)) be initial conditions satisfying 

(4.18) sup(re h + \u h (t )\ + K(*o)| + \u h (t )\ + |f h (*o)|) < +oo. 

h 

Then the discrete d'Alembert-Lagrange-type principle (I4.14p relative to Oh has a 
unique solution (uh,Vh)- 

Moreover, there exists (u,v) G W 2, °°(I, M. DxM ) satisfying (|4.8j) such that, up to 
a sub-sequence, ((uh, Vh))heN converges to (u, v) uniformly on I , and ({iih, f h))/ieN 
converges to (u,v) in L p (I,M. DxM ) for every p G [l,+oo). 

To prove the theorem above, preliminary results are necessary. As a first step by 
Lemma [T] we establish an explicit form of the discrete variational principle (|4. 14[) . 
From such a form, it is easy to infer existence and uniqueness of discrete solutions 
with given initial conditions. Then we show BV estimates for the velocities of 
discrete solutions in Proposition [1] 

By making variations that fix all the nodes except one at a given time, one can 
find easily the following equivalent form of the d'Alembert-Lagrange-type principle 

Lemma 1. Let (u,v) G y@, then (u,v) solves (|4.14[) if and only if given any node 
a G T and any nodal time t\ G (to,tf], denoting by K the sole element for which 
t l a G Ok and t K = t l a , if a G T\ {dB u U dB v ), the following balance equations hold: 

m a (u a (t*- 1 ) -u a (4)) 

= (*K 1 - 4) / d F e ((VU K , V K , VV K )(X, 4)) VAf a (x) < 



I dx 



J d u w([u K ,u K ){x,i J K ))M a {x)dx 



(4.19) t{x)N a {x)dU d -\ 



dKndBt 



an 



d 



Pa (fa (C 1 ) - »« (4)) + *7>a(4 +1 ) - f .(**)) 

= ^k 1 -^k)J i d N e({Vu K ,u K ,Wv K ){x.t ] K )^Ma{x)dx 

(4.20) + (i J A + 1 - 4) / d v w f*)C=>**)) AA q (x) ■ 



i dx. 

Moreover, if a G dB u \dB,j, only the first equation has to be satisfied, while, if 
a G dB v \dB u , only the second equation has to be satisfied. 
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Clearly, equations (|4.19|) and (14.20P define inductively a unique discrete trajec- 
tory in 3^e once initial conditions have been specified. 

Estimates on the L°° norm on the time interval I of the velocity of stationary 
points can be derived by exploiting directly the discrete stationarity conditions 
and the growth conditions of the potential energy densities as suggested in 
(see also [8]). Here, we provide also an estimate on the pointwise variation of the 
velocities of stationary points which has been overlooked in previous works. In 
turn this estimate let us gain strong compactness properties and permits us to 
pass to the limit directly into equations (|4.19p and (|4.20[) , or equivalently in (|4. 14[) . 
This method is alternative to the ones already developed to study the convergence 
properties of variational integrators (see [TO], [19], [11], [8]). Furthermore, the BV 
estimate wil be instrumental in Section [5] (see the related comments there). 

Proposition 1. There exist a constant Tq — Tq(t], p, A, A, T) > such that given 
initial conditions (u(io), ^(^o)) an d {u{to),v(to)) for every entire time set with 
Tq G (0,To), the solution (u, v) 6 J^e of equations (|4.19|) and (|4.20| satisfies 

(4.21) IHIl~(/,r«) + IHU°°(/,rm) < «exp (kt&) , 

for some constant n > depending on the initial conditions themselves and on the 
data of the problem. 

Moreover, the functions u, v have (pointwise) bounded variation pV on I , with 

(4.22) pV(u, [ti,t 2 ]) + pV{v, [ti,h]) < KT e C 1 + exp (zero)) (*a - *i + 2T e ), 
for every interval C J. 

Proof. Fix a node a in the tessellation T and a nodal time t* 6 Q a . By assumption 
there exists a unique K S T such that t l a £ Qk, with t l a — t J K . The definition of 
nodal time set yields t l a = ta"^ G Q a for all the other nodes a e K. Moreover, 
recalling that 6 = {i/};=i,...,Ar e , we assume t l a — t s+ i e 9 for some s. 

We first estimate the velocity of u at time t s+ i in terms of the velocities of u itself 
and v at previous times. In doing that we use the discrete d'Alembert-Lagrange- 
type equation (|4.19[) , assumptions (A1)-(A3), and (|4. 6[) which entail 

(4.23) m a \u a (f a )-u a (C 1 ) I < c(t^-t^ K ) ( 1 + £ (\u a (t^)\ + \u a (^)\)) . 
By taking into account that the maps u a , v a are piecewise affine in time, we get 

\u a (t*)-« a (c 1 )! <c(4 +1 -4) ( 1+ E (K(*o)i + K(*o)i)j 

E (4 +1 -4)(M4)l + M*L)l) 

< cT e + cT e (t f -t ) {j^h^^ta.tT 1 - 1 )^) + "^^-((to.tL^- 1 )^)) ■ 

By recalling that = t s +i, and since ia < is for every a 6 Jf, the latter 

inequality yields 

(4.24) 

INL»((w«+i),K D ) - cTe ( x + W^h^ato^)^)) + (i + cT e) l|w|lio O((40ita)iRD) . 
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We prove the analogous estimate for z>, the only difference with the previous argu- 
ment being that we have to take into account, in addition to the previous case, for 
the presence of the discretized dissipation. Thus from (|4.20p . assumptions (Al)- 
(A3) and (|4~6|) we get 
(4.25) 

~Pa\K (4) - K (C 1 ) I < VMu a (4) I +cT e ( 1 + (K(*a (0) )l + I "«(*a ( °° ) I) J . 



aeK 



which in turn yields 

k (4) - K (C 1 ) I < cT e \v a (4) | + cT e 



+cTe(t f -t ) 



Thus, there exists a constant To = Xb(?7, p, A, A, T) > such that for every Tq G 

(0, T ) we have (1 - cT e ) > 1/2. Finally, we deduce 

(4.26) 

IWW L^((t (h t B + 1 )M M ) ^ cT © i 1 + IHlL~((t ,t s ),K D )) + ( X + cTe ) ll^lli^CCio.t.J.R-M) • 

In pa rticular, by setting /?, := ||u||i«,(( tDjtl))R r>) + \W\\L^atoM)M M )^ from 1 4 - 24 ^ 
(14361) we infer 

/? fl +i <cr e + (i + cr e )/3 s - 

Such inequality yields by iteration 

s 

(4.27) /3 S+1 < cT e £ (1 + cT e ) 1 + /?„ (1 + cT e ) s+1 

t=o 

< (1 + p ) (1 + cT e ) s+1 <c(l + cT e f e . 

Eventually, since Tq#Q < (if— to)T&jfcT, with t@ defined in (|4.1ip . from inequality 
P~27) we deduce 

||w|U«>(j,m D ) + IHIl~(/,r m ) =/9#e < c (1 + cT e ) cre/T ° , 

and then (|4.21j) follows. 

In order to prove the BV estimate (I4.22p we note first that since ii a , v a are 
piecewise constant, their (pointwise) variation over an interval [ii,i2] C J is given 
by the sum of the jumps in [ii,t2]. Hence, by using (|4.23|) and (I4.25| . and taking 
advantage of (14.211) . we obtain 

(4.28) \u a (4) - ii a (C 1 )! + \i>a (4) - K (C 1 ) | < cT e (1 + c exp( C T e ) ) . 

To compute the pointwise variation of u a , v a on [ti,<2] we sum (|4.28l) over the set 
of indices S tljt2 := {i\ [t' l a ~ 1 ,t 1 a ] n [^1,^2] 7^ 0}- Furthermore, by taking into account 
that Te#S tli t 2 < ere (£2 — ii + 2Te), we conclude that 

= E (M^-Mc^l + M^-Mc 1 )!) 

ies tl , t . 2 

< cTe#S tl)t3 (l + cexp(cre)) < kt b (1 + Ke^p(KT e ))(t 2 - ii + 2T e ). 

□ 
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Lemma 2. Let (@h)heN be a sequence of time sets such that Te h — > as h — » +oo. 
Let {{vhi Zh))hen C M DxM ) converge strongly to some (v,z), then 

(4.29) \im At, e h (vh,z h ) = A T (v,z), 

h 

and for every (ip h ,(f) h ) G ^e h H Hq(I,M. DxM ) swc/i i/iai ((V/t) converges 
strongly to some (ip,4>) G i? 1 (/,R D><M ) 

(4.30) lim<5^r,o fe (w/ l ,2/ 1 )b/ l ,0/ l ] =MT(«,2)[ft^]. 
Moreover, it holds 

(4.31) lim2? T) e,(4)[^]=2?r(i)M. 

Proof. The convergence property stated in (|4.29|) follows easily from the strong 
convergence of ((wh, -z/i))/ieN iri H 1 and the regularity conditions of the integral 
densities of V assumed in (A1)-(A3). 

Define ip h := (v h ,z h ), ip := (v,z), $ h = ((p h ,<f> h ) , $ : = {(f,4>). Let (Sj) be a 
positive vanishing sequence. By taking into account (|4.29|) . for every j G N we get 

h \ Sj Sj J 

Hence, with fixed any subsequence (hi), by a diagonal argument we find a further 

subsequence (hi . ) such that 

(4.32) 

' ^r,@ h ..{ip hi . +5^ hi .)-A r ,& h .i^ hi .) A T (ip + - A T (ip)\ n 



lim 

3 



Notice that by definition of first variation 
(4.33) 5AT W m=^^ L+m ' ATW 



j Si 



= 0. 



J 

Moreover, for every h the action Ar,Q h is a C 1 functional over y@ h by the regularity 
assumptions (Al)-(A3) on the densities of the potential V. Thus, for some Sj G 
(0, Sj) we have 

Ar,e hi . (i> hi . + Sj$ h ) - A T ,e h . . (ip hi . ) 
(4-34) '- '- -=-? ^ 2- - Mr,8 k , W> hi . 

Oj 3 3 

Actually, the same regularity assumptions imply that 
(4.35) lim(sA T ,e, H .(^ hz . + e^hj^hj - SA r ,e h ._ (Vv ]) 

J \ 3 3 J 3 3 J 3 3 ) 

Thus, by (|4.32p - (|4.35p we infer (|4.30p for the subsequence (hi,). Urysohn property 
justifies (|4.30[) for the whole sequence. 

To prove (|4.31[) , note that the semi-discretized dissipation XV is continuous along 
sequences strongly converging in H , so that 

limV T (z h )[<P h ]='D T (z)[ct>]. 

h 

To obtain f)4.31j) we need only to prove 

lun(V T ,e h (zh)lM ~ V T (z h ){M) = 0. 

h 

The limit above is a consequence of (14. 17)) (since Te h — » as h — > +oo) and the 
strong convergences in H 1 of (zh)he® and (4> h )h&{ to z and 4>, respectively. □ 
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Proposition Q] and Lemma [2] are necessary tools for the proof of Theorem [2j 
Proof of Theorem [H According to (|4.18p we may fix t > such that for every h 

T~e h + \uh(to)\ + \vh(to)\ + \uh(t Q )\ + \vh(to)\ < tq. 

Denote by (u^, vh) the solution to equations (|4.19|) and (|4.20[) . associated with O/j. 

By Proposition [TJ we infer that ((uf l ,i'h))heN is bounded in BV(I,M. DxM ) and 
by the BV compactness and embedding theorems (see [U Theorem 3.23], p] Corol- 
lary 3.49]), we deduce the existence of functions (u, v) G W 1 '°°(I,R DxM ) such 
that (up to a sub-sequence) ((uh, Vh))h£N converges to (u, v) uniformly on I, and 
((iih, Vh))h£H converges to {ii,v) in L p (I,R DxM ) for every p G [1,+oc). 

Actually (u,v) G W 2 >°°(I, R DxM ) as observed for the one-dimensional case in 
[11)1 Theorem 4.7]. We repeat here the proof for the sake of completeness. According 
to (14.28[) . for every node a G T and for every t\, t 2 G I we have 

\{u h ) a (h) - (U h )a(tl)\ + \{Vh)a{h) - {v h )a{tl)\ 

(4.36) < ktq (1 + exp («r )) (t 2 -h + 2T e J. 

Fix a standard mollificr tp G Cq°(— 1,1), with ijj > and f.i/jdx = 1 and set 
ip s (x) = 5~ 1 ip(x/5). Let (u s h , v s h ) — (u h * tp g , v h * ip s ), it follows from (|4.36D 

\(u h ) 5 a (h) ~ (u h ) 5 a (h)\ + \(H) 5 a (t 2 ) - (p h ) s a (h)\ 

< ktq (1 + exp(KTo)) (t 2 - ti + 2T e J + cS. 

By taking into account Lebesgue point theorem for (u a , v a ) and passing first to the 
limit as h — > +oo (note that by assumption Te h — > as h — > +oo) and then as 
(5 — > in the inequality above, we conclude that 

\u a {h) - u a {h)\ + \v a (h) - i>o(*i)| < kto(1 + exp(KT ))(£ 2 - h). 

Eventually, we show that (u, v) solves (|4.8[) by passing to the limit in (14. 14j) . Take 
test functions (<p,</>) G H^(l,R DxM ), and let (fp h ,4> h ) G y@ h H H^(I,R DxM ) 
be a sequence of piecewise affine test functions converging to (ip, <fi) strongly in 
H 1 (/, R DxM ) (see for instance [THl Lemma 4.3]). Then the discretized d'Alembert- 
Lagrange-type principle (|4. 14)) is satisfied, namely 

^r,e h (u h ,v h )[<p h ,^ h ]+VT,e h (vh)[i> h ] = 0. 

The conclusion then follows by Lemma [2] and the strong convergence in H 1 of all 
the relevant sequences. □ 

Remark 1. A standard diagonalization argument implies that we can allow for 
intervals I = (to,+oo) provided that all the statements in Theorem^ are intended 
in a local sense. 

Remark 2. The methods developed to prove Theorem^ can be refined to deal with 
kinetic co-energies \ of the form 

X{ v ,v) = \^ ■ " 

where Q is a positive definite symmetric element of Hom(T u M.,T*M.). The idea is 
to choose for the fields u, v different shape functions: for u the family {Af a }aeT as 
above, and for v a family {A/"„}aeT in such a way that the mass matrix associated 
with the spatial discretization of \ is in diagonal form. 
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5. Synchronous Variational Integrators 

Although the choices \pv° ■ v and \£lv ■ v as expressions of the substructural 
kinetic energy cover a wide set of special cases, there are circumstances in which 
more complicated expressions arise (see for example the dynamics of liquids with 
a dense distribution of bubbles [5]). For this reason the analysis of the dicrete 
schemes in which the substructural kinetic co-energy is left unspecified has its own 
importance. In this case we prove the convergence of the discrete dynamics only 
for synchronous variational integrators, i.e. when @k coincide with © for every 
K e T. 

More precisely, we assume that x is anv function with the following properties 

(A4) x € C 2 (M. kxk , [0, +oo)) has quadratic growth, i.e. there exists a constant 
S > such that for all v, C € R k 

|v 2 x(^,C)l<s. 

(A5) x(y-i ') i s uniformly convex with respect to v, i.e. there exists a constant 
7 > such that for all u, (, £ E R k 

The spatial discretization follows the same lines indicated in Subsection 14.11 In 
the sequel we use the same notations taken there, the only difference being that 
the family of shape functions {Af a } a &T is chosen to form an orthonormal system in 
L 2 (B) (see Remark[3]). A discrete version in space of the d'Alembert-Lagrange-type 
principle (|4.ip arises. We then get 



m a u a (t)-<p a (t)dt = J d F e((Vu,v,Vv)(x,t))VAf a (x)dxj-<p a (t)dt 
(5.1) +1(1 d u w{{u,v){x,t))M a {x)dx- f t(x)Af a (x)dH d - 1 (x) ) -<p a (t)dt, 



and 



di,x((v,u)(x,t))J\f a (x)dxj -<f> a (t)dt 

d v x((v,i>)(x,t))Af a (x)dx) ■ 4> a {t) dt + r, [ v a {t)-<t> a {t)dt 
: / J i 

d N e ((Vu, v, Vv){x, t)) VAfa(x) dx ) ■ <j> a (t) dt 
(5.2) + I ( [ {d v e{{Vu,v,Vv){x,t))+d v w{{u,v){x,t)))N a {x)dx\ -(p a (t)dt, 



where the test functions satisfy (p a E HQ(I,M. d ) and tp a = if a E dB u , and 
(f> a E H$(I, R k ) and (f> a = if a E dB u . 

As regards the time discretization, we fix a discrete time set = {ii}i=i ) ...,jv e 
for /, with time size ft-e — maxe(ti+i — tj), and consider the subspace 3^e of 
H 1 (I,M. DxM ) of continuous functions which are affine on each interval (£j,ti+i) 
and such that u a = if a G dB u and v a = if a £ dB u . We then define the discrete 
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action as 

i=0 

+oc if ey\y e , 



At,® (u,v) 
where 



A l Te {u,v) := V \\u a {U)\ 2 + / x((i/,*)(z,t i ))£*c-V'((u,i')(ti))i 
aG r 2 ^ 

and the potential V is defined in (|4.3[) . Notice that in this case 0if = for every 
K £ T, and thus To = ^e, T e — 1- We denote by r e the ratio 

( 5 - 3 ) r e : = ~n 77 7T- 

mine {U+i - U) 

Note that we have chosen the approximation 

rU+i / r \ 

X{{v,v){x,t))dx-V{{u,v){t)))dt 



In this setting, the d'Alembert-Lagrange-type discrete principle is then given by 

(5.4) Mr,e(u,v)[M]+2V,e(£')M = 0, 

where (93, 0) £ 34s n if 1 (/,R £,xM ) and £> r . e is defined as in (j4~T5]) . namely 

5r,eWM == »» S 53K(*0 - "«(<i-i)) ■ 

oeT z=l 

The main result of the section reads as follows. 

Theorem 3. Assume (A1)-(A5). Let (Qh)heN be a sequence of time sets for a 
bounded time interval I — (to,tf) such that Tg> h — > as h — > +00. Let also 
(uh{to), Vh{to)) and {uh{to) , v h{to)) be initial conditions satisfying 

(5.5) sup (r' eh + \u h (t )\ + \v h (t )\ + \u h (t )\ + \vh(to)\) < +°°- 

h 

The discrete d'Alembert-Lagrange-type principle (|5.4|) relative to Oh has a unique 
solution [uh, Vh) for h sufficiently large. 

Moreover, there exists (u,u) £ W 2 '°°(I, R DxM ) satisfying (|5.ip and (|5.2p such 
that, up to a sub-sequence, ((«/,, Vh))h&1 converges to (u,i>) uniformly on L, and 
((uh,i'h))heN converges to {u,v) in L p (I,R DxM ) for every p £ [l,+oo). 

The proof follows the same strategy developed for the proof of Theorem [5J once 
all the preparatory steps have been accomplished. Thus we will not give the details 
of the proof of Thcorcm[3l Instead, we will focus our attention to show existence and 
uniqueness of discrete solutions to (|5.4|) and the BV estimates on their velocities. 
In doing that the main ingredient that we will exploit is the uniform convexity of 
•) assumed in (A5). 

In the following lemma we collect all the properties of x which will be useful in 
the sequel. The proof of the lemma is standard and then omitted. 

Lemma 3. Assume (A4) and (A5), then the following properties hold true 



DISCRETE DYNAMICS OF COMPLEX BODIES 



19 



(1) Uniform Strict Monotonicity of di,x{ v > •)• F° r a ^ u > Ci> C 2 

7IC1 - C2I 2 < (d>x(y, Ci) - hxiy, C 2 ). Ci - C 2 >- 

(2) Uniform Lipschitz Continuity of <9,>y(-, f^xC^i ■)•' ^ or a ^ ^2, C G ^ fc 

l^x(^i,C) - #<>x(^2,C)l < S K - 

and for all v, Q 1 , ( 2 G K & 

l^xKCi)-Svx(f,C 2 )l<s|Ci-C 2 |. 

(3) Growth conditions of Xi ®vXi &vX : There exists a constant S2 > such 
that for all v , £ € K fe 

|^x(i/, C)| + |^X(^C)I <H 3 (1 + M + |C|), 

7 |CI 2 - s 2 (i + |H 2 ) < x(^,C) < Sa(l + M 2 + ICI 2 )- 

5.1. Convergence in time. As a first step we write an equivalent form of (|5.4|) 
which can be obtained by varying only one node at each time G O. 

Lemma 4. Let (u,v) S ^e- solves (|5.4p i/ and oniy i/ given any node 

a e T\(dB u U 9,6„) and any tone £j £ H (io, i/], £/ie following balance equations 
hold: 

m a (u a (ti-i) - ii a (U)) = (t i+ i - U) / d F e{{Vu,v,Vv){x 1 t i ))VJ\f a (x)dx 

JB 

(5.6) +(t i+1 -ti)( [ d u w((u,v)(x,ti))N a (x)dx- [ t(x)Af a (x) dH d ~ x 

\JB JdBt 

and 

{dbXiivityfaU-i)) ~ dox{{v,v){x,ti))) N a {x) dx 



+ (ti+l~ti) / d v X {{V, v){x, U))N a {x) dx + T)(v a (ti+l) ~ Va(U)) 
JB 

= (t*+i — *i) / d N e((Vu,v,Vi'){x,t i )) VN a {x)dx 

JB 

(5.7) +(ti+i-*0 / (ft,e((V«,v,Vi/)(ar J * < )) + ft / to((«,i/)(aj,t i )))A/'„(a:)d!B. 



In case a € dB u \dB u only equation (|5.6p /las to fee satisfied, while if a 6 dB u \dB u 
the sole (|5.7p /ias to fee satisfied. 

If for a fixed time tj 6 On [to>*/) a solution to (|5.6p and (|5.7p exists for every 
node in T, a unique solution at time ti+± is determined. An induction argument 
provides the conclusion once the initial conditions are specified. 

Consider the function * : Rkx#(T\dB„) R fcx#(T\3B„) ; ^ = (tf a ) oeTV , Bi/ , 
where for every node a £ T \ 9£>„ 



*«M := / d^x\ ^{x,ti),J2M a ( 



— V<xt% ) Af a (x)dx 

H+l — H I 



d i ,x{{v 1 v)(x 1 t i ^ 1 ))Ma{x)dx + n{v a - i/ a (*i)) 



(5.8) +(U+i-U) I d vX ( v(x,U), y ] M a {x) — Va{ \ l) \M a {x)dx. 

H+l — H J 
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Notice that for every node a S T \ dB v equation ()5.7j) becomes 

* a (("a(*i-t-i))aeT\aB„) = (U+x ~U) / 9ive ((Vit, i>, Vi/)(ar,tj)) VjV a (x)dx 
(5.9) - tO / (ft,e((V«,i/,Vi/)(a! > * < )) + ft,to((u,i/)(a:,t i )))Ar o (x)da:. 



Recall that v a (t) = for every time t S 7 if a € <9£>i,. 

Lemma 5. There exists T — Tq(7, 3) > suc/i t/iat, /or any time set G with 
Tq € (0, To), t/ie map 5" is a bisection. Moreover, for every v,fj, £ R fcx #( T \ aB ") 

(5.10) 2 ^ |i/ tt -/i a | 2 <<*W-*(M)^-A*>. 

Proof. Let us show the injectivity of VE". Given ^ and \x in R fex #( r \ 98 ") ; the uni- 
form strict monotonicity of di,x(v, •) (see (1) in Lemma [3]), the uniform Lipschitz 
continuity of d^xi'^) ( see (2) in Lemma [3]), and the orthonormality of {Af a }aeT 
in L <(B) imply 

(V(v) - *(/*), V - fi) = <*aM-*aO*).fa- Ma) 

' 77 |l/a -^ |2+ (^f-^- s ) E K "^' 2 



> 



If To < 27/(2S + 7), the latter equation implies inequality (f57TO|> for T e £ (0,T ). 
The injectivity of ^ is then established. 

Eventually, the surjectivity of W follows at once from (|5. 10|) by arguing compo- 
nentwise. □ 

Existence and uniqueness of solutions to (|5.6p and (|5.7p can be shown. 

Proposition 2. //To = Xb(7,S) > is the constant in Lemma\5\ for every time 
set with size Tq £ (0,Tb), a unique solution to equations (|5.6| and (|5.7p exists 
in y® once initial conditions (u(to),v(to)), (u(to),P(to)) are given. 

Proof. By taking into account (|5.9p and the invertibility of <J/ proved in Lemma [5] 
(for Te £ (0, To)), (I5.6[) and (|5.7p define inductively a unique piecewise affine tra- 
jectory (u, v) G y@ once the initial conditions are specified. □ 

It is clear that in this setting the variational approach based on T-convergence 
proposed in [TH] cannot be pursued. In the conservative case a more conventional 
approach based on the repeated use of Gronwall inequality has been developed in 
[TDJ. 

The following proposition is the analogue of Proposition [TJ Let us point out that 
the BV estimates proven in the sequel are instrumental. Indeed, they enable us to 
pass to the limit directly into the discrete equations (|5.6p and (|5.7p and to prove 
that their solutions have cluster points solving (|5.ip and (|5.2p . 
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Proposition 3. There exists a constant T\ = 7i(7, S, r), A, A, T) > such that 
given initial conditions (u(to), v(to)) and (ii(to), i'(io)) for every entire time set 
with Tq £ (0,Ti), the solution (u,ii) £ y@ to equations (|5.6p and (|5.7p satisfies 

(5.11) ||«IU~(7,kd) + IHU°°(/,rm) < Kexp(KTe), 

/or some constant k > depending on the data of the problem and on the initial 
conditions themselves. 

Moreover, the functions it, v have (pointwise) bounded variation on I with 

(5.12) P V(u, [t 1} t 2 ]) + P F(*>, [ti, t 2 ]) < kt'q (1 + exp («r' e )) (i 2 - f x + 2T e ), 
/or every interval [tijia] C J. 

Proof. Fix a time i; G ©. To estimate the velocity of u at time U in terms of 
the velocities of u itself and v at previous times we follows the path indicated in 
Proposition [TJ We use equation (|5.6j) and, for every node a £ T, we get 

(5.13) |u (^) - w Q < cT e (1 + ||u||i«»((t ,t < _ 1 ),K^') + IHU~((t ,t*-i),KM)) ■ 

The latter inequality yields 
(5.14) 

lL~((to,ti),R D ) - cTe {} + H^lli^^to^j-i),^)) + (1 + cre ) MIl°°((W;-i),R d ) • 



To prove the analogous estimate for v we use the equivalent form of equation ([5 
written in (|5.9[) . We rewrite the left hand side in (|5.9p as 

*a((^a(^+i))aer\aej = + 4 2 + ^a, 

where 

J a := / (^XK 1 !*!)!^)*!-!)) ~ d i ,x{{v,v){x 1 t i ^ 1 )))M a (x)dx, 

Ja '•= -(*<+!—*») / ((^^OM*)) A/" a (l) cfe + 7?(^a(*i+l) - ^a(U)). 



To evaluate the variation of z> we consider the scalar product of \& calculated in 
(^a(*i+i))aer\aB„ with the vector 77 := (i/ a (ti) - v a (U-i)) ae T\dB v > so tnat we get 

(*((^ ct (^ + i)) Qe r\a8j,f 7 ) = ^{Jl^aik) - v a (ti-{)). 

r=l,2,3 aeT 

Indeed, recall that v a {ti) = for every i if a £ <9£v. We estimate separately each 
term above. First by (1) in Lemma [3] we have 

-*a(«i-l)> = 

(^x((^»(x,ti)) - (^(x,ti),P(x,ti-i)) , v(x,U) - ^(xjti-i)) dx 
(5.15^ 7/ |z>(o;,ij) - v{x, £j_i)| 2 dx = 7 V] |/> a (*<) - ^a(*i-i)| 2 = l\v\ 2 ■ 
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Then, by (2) in Lemma [3] and since u(x,ti) — u(x,ti-i) = v(ti-\){ti — tj_i), by 
Holder inequality we obtain 

= / {duX ( v ( x , U), v(x, U-i)) - dvX ((v, v)(x, U-i)) , v(x, U) - v(x, ij-i)) dx 

JB 

(5.16) -3T e ||£'|| £ ,oo(( i0)t4 _ 1 ) ) ] R M)|F|. 
Moreover, by (3) in Lemma [3] we infer 

5Z(^a; Va(U) - U a {U-l)) = 1 - U) ^(^0(^)1 ^a(k) - V a (ti-l)) 

aeT a£T 

-(ti+i-k) / {d v x{{v,v){x,ti)) ,v{x,ti) - u(x,ti-i)) dx 

JB 

(5.17) > -cT e (v + =1) (l + IMlL~((t ,*o,R M ) + IHU°°((i ,t«),R M )) \ v \- 

Furthermore, we use (A2) and the fact that d v e, d^e have linear growth to bound 
the right hand side of (|5.9p as follows 

(5.18) < cT e (1 + ||u|Uoo ((t0)ti _ l)iR D) + |HU=o ((Wi _ l)iKM) ) \V\. 
By collecting (f3TT5]) - (f5TT5|) it follows that 

(5.19) 7|F| < cT (l + ||i'|U»((t 0l t 4 ),K-» f ) + [l li llL«((to 1 ti-i),B D )) . 
which in turn yields 

(7 - c7 e)||t'IU»o((t ,t < ),RM) 
< c7l|i / lli»((t ,* i -i),R M ) + cT e (1 + |[«IU«'((to,ti-i),m D )) • 
Thus, there exists a constant Ti = ^(7, S, 77, A, A, T) > such that for T e £ (0, Ti) 

(5.20) IHIiooj^.^RM) < c||z>|| L = ((t0iti _ l) ^M ) + cT e (l + Wuh^UtaA-i),® )) ■ 

In pa rticular, by setting j3 l := ||u||i«,( (t0jtl))R r>) + [[^Hz-o^^rm), from (|5.14[) and 
(15T2T))) we infer 

ft <cT e + c(l + Te)^_ v 

By iterating this inequality and arguing as in the proof of Proposition [I] one can 
deduce the L°° estimate (15. lip . 

Eventually, in order to prove the BV estimate (I5.12j) we can reason as in the 
proof of Proposition Q] by taking into account (|5.1ip , (15.131) and (|5.19[) . □ 

Remark 3. The orthogonality condition of the family of shape functions {A/" a }oeT 
has been imposed only for convenience. Indeed, the same model of Section [^J in 
which the potential is decomposed element-by-element, can be developed for the syn- 
chronous setting, too. 
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